function F = PredIntegral1(statec, absBank, StogoExp, tgoExp)
global Vf N dtInt alphaLD ef Vup Vdown alphamax
rad = statec(1);
V = statec(4);
fpa = statec(5);
dtIntPred = tgoExp/N;
if dtIntPred <= dtInt
    dtIntPred = dtInt;
end
stateLon = [rad V fpa StogoExp];
while 1
    V = stateLon(2);
    if V > Vup
        alphar = alphamax;
    elseif V > Vdown
        alphar = (alphamax - alphaLD)*(V - Vdown)/(Vup - Vdown) + alphaLD;
    else
        alphar = alphaLD;
    end
    ctrl = [alphar absBank];
    stateLon = RK4(@LonPredDynamic, stateLon, ctrl, dtIntPred);
    e = 1/stateLon(1) - stateLon(2)^2/2;
    if e >= ef % stateLon(2) <= Vf
        break;
    end
end
F = stateLon(4);
end